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We explore a scenario of New Physics entering the description of B —> K 1 '*' 1 pp decay through 
couplings to the operators Og^o, satisfying Cg = —C[ 0 . From the current data on B(B 3 —> pp) 
and B(B — > Kpp) ns^lGeV 2 > we obtain constraints on ReC[ 0 and ImC[ 0 which we then assume 
to be lepton specific, and find Rk = B(B —> K pp)/B(B —¥ Kee)^ ,e]G eV 2 = 0.88(8), consistent 
with recent value measured at LHCb. A specific realization of this scenario is the one with a scalar 
leptoquark state A, in which C[ 0 is related to the mass of A and its Yukawa couplings. We then 
show that this scenario does not make any significant impact on B s — B s mixing amplitude nor to 
B(B — > Kuv). Instead, it can modify Rk * = B(B — > K*pp)/B{B —> A*ee)n 6 ] GeV 2 , which will 
soon be experimentally measured and we find it to be Rk * = 1.11(8), while Rk* /Rk = 1.27(19). 

A similar ratio of forward-backward asymmetries also becomes lower than in the Standard Model. 

PACS numbers: 13.20.He,12.60.-i,14.80.Sv 


I. INTRODUCTION 

The b —> s transitions were in the focus of many theoretical and experimental studies during the last two decades 
due to the possibility to constrain potential New Physics (NP) contributions at low energies. With LHC7 and LHC8 
runs direct searches for NP became available. This gives us an excellent opportunity to question the appearance of 
physics beyond Standard Model (SM). At low energies B- factories and the LHCb experiment provided flavor physics 
community with a lot of rather precise results on b —> s transitions. The LHCb experiment has observed slight 
discrepancies between the SM predictions and the experimental results for the angular observables in B —> K* p + p~ 
decay. This effect has been attributed to NP, although the tension might be a result of the SM QCD effects. Recently, 
another anomaly in b —> s£ + C~ transition has been found in the ratio of the branching fractions, 

_ B(B — > Kp + p~) q 2 e [ lfi ] GeV 2 

B-K ~ m r> - , _ N - ■ l 1 ) 

B\B -)> Ke+e ) g 2 6 [ lj6 ] GeV 2 

LHCb Collaboration measured this ratio for the square of dilepton invariant mass in the bin 1 GeV 2 < q 2 < 6 GeV 2 , 
and found [1], 

i?LHCb = q 745 ± 0 ; 090 -1-0.036 , (2) 

lower than the SM prediction, RfP 1 = 1.0003 ± 0.0001, in which next-to-next-to-leading QCD corrections have been 
included [2]. In other words, the LHCb result points towards a 2.6 a effect of the lepton flavor universality violation. 

Furthermore, the combined data analysis of the B s —> p + p~ events gathered at LHCb and CMS resulted in B(B S —i 
p + p~) = (2.8lo'g) x 10~ 9 [3]> in good agreement with the SM prediction B(B S —i p + p~) = (3.65 ± 0.23) x 10 -9 [4]. 
This offers an excellent probe of b —> sp + p~ transition in the light of SM and gives rather tight constraints on 
parameter space of many models of NP. The Rk anomaly has been approached in the literature in different ways: 
either by using the effective Lagrangian approach or in a specific model of NP. For example the effective Lagrangian 
approach used in references [5-8] indicated that in order to understand the measured value of Rk one must include 
the effects of NP, and that the effects of non-perturbative QCD alone could not explain such a large deviation of Rk 
from unity [6-21]. In particular, it was found that the NP contribution most likely affects Cg, C±o or C' g , C' w effective 
Wilson coefficients, and that some kind of lepton flavor universality violation is needed, e.g. C g ^ Cg [10, 22]. In 
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order to determine whether Rx anomaly is due to NP in electron or/and muon couplings through a combined analysis 
of several decay modes, it is very important to have a high precision knowledge of hadronic form factors [16-18], which 
can be computed in the region of large q 2 ’s by means of numerical simulations of QCD on the lattice [23-25]. 

In this study we first use a model independent approach, assuming that NP contributes at low energies to an 
operator that is a product of a right-handed quark and a left-handed muon current. In the language of b —> s/x/x 
effective Hamiltonian such a situation corresponds to a combination of Wilson coefficients Cg and C' w , and that they 
obey Cg = —C[ 0 . Decays to the final states with electron-positron pair are instead governed by the SM only. This 
assumption is motivated by the fact that measured quantities of b —» se + e~ processes agree with the SM predictions 
better than they do for the b —» s/x + /x _ processes [12], which are also more precisely measured than the electronic 
modes. We consider simultaneously the constraints posed by B(B —» AT/x + /x“) and B(B S —» /x + /x“) on such a scenario, 
and then predict the Rx as well as Rx * ■ We discuss other observables which might serve as additional probes of the 
observed lepton-flavor universality violation. 

A specific realization of the scenario we discuss in this paper is a model with a light scalar leptoquark A with 
quantum numbers of SU(3) C x SU(2)l x U(1)y being (3,2,1/6). It indeed verifies the relation, C' g = — C' w [9], 
and leads to a consistency with the measured value of Rx- The features of this leptoquark state have been already 
described in the literature [26] . While there is no theoretical motivation to forbid leptoquark contributing to b — > see 
decays, simultaneous presence of both muonic and electronic couplings could be problematic because they would, 
together, induce lepton flavor violation in B s —> eg, and /x —> ey decays. It is interesting that the flavor physics 
constraints at low energies agree and are complementary with the constraints obtained from the direct experimental 
searches at LHC [27, 28]. Furthermore, the atomic parity violation experiments provided a strong constraint on the 
interaction of the down-quark-electron interaction with the leptoquark state [26, 29], while the couplings to muons 
appear to be less constrained via B(K R ► g ± e zf ) < 4.7 x 10~ 12 [26, 30]. We therefore assume in our analysis that 
in the b —> st + £~ processes only the muons can interact with the leptoquark state. A few other leptoquark states 
have been discussed in the literature [6, 9, 14, 16] as possible candidates to contribute to the Rx anomaly. However, 
the leptoquark with quantum numbers (3, 2,1/6) has a desired feature that it can be light without destabilizing the 
proton [31-33]. Notice also that another light leptoquark scalar state, not mediating the proton decay, is (3, 2,7/6) 
and it leads to the relation Cg = C\g. That latter scenario, however, cannot explain the Rx anomaly as discussed in 
Refs. [6, 14]. 

In Sec. II we remind the reader of the main definitions and give basic expressions for B(B S —> g + g~) and B(B —> 
Kg + g~), which are then used, together with the experimental data in Sec. Ill, to constraint C' l0 = —Cg and show the 
consistency of our value for Rx with the measured one at LHCb. Furthermore, we make a prediction of the similar 
ratio in the case of B —» K*l + £~ decays and discuss other observables that might be of interest for testing the lepton 
flavor universality violation. In Sec. IV we discuss a model with scalar leptoquark in which the relation C' w = —C'g 
holds exactly, and is connected to other similar processes involving the b s transitions which we also discuss. We 
finally summarize our findings in Sec. V. 


II. EFFECTIVE HAMILTONIAN AND BASIC FORMULAS 


The processes with flavor structure ( sb ) (/ j,g ) at scale g = g b = 4.8 GeV are governed by dimension-6 effective 
Hamiltonian [34-36] : 


B-cS 


4 G f 


v tb v* 


6 

Y.CiMOifa) 

1=1 


E (Ci^o^ + cK^o'M) 

i— 7,...,10 


( 3 ) 


The contributions of the charged-current operators 0 1 , 2 , QCD penguins © 3 ,. ..,6, and the electromagnetic (chromomag- 
netic) dipole operators O 7 (Os) will be assumed to be saturated by the SM. On the other hand, operators involving 
a quark and a lepton current will contain the SM and potential NP contributions. The basis of operators may be 
further extended to account for possible (pseudo)scalar or tensor operators [23], whereas for the purposes of this work 
the following operators will suffice: 

O7 = Arm b (s(J llu P R b)F> 1 '' , C > 8 = -m b (sa^G l "'PRb) , 

9 2 9 2 ( 4 ) 

0 9 = ^(ry ll P L b)(i'fe) , 0i„ = • 

g 2 5 

Here Pl/r = (1 =F 7s)/ 2, while e is the electromagnetic and g the color gauge coupling. F and C 1 ' are the 
electromagnetic and color field strength tensors, respectively. The basis is further extended by the wrong-chirality 
operators, O'g^g, which are related to Og^g by replacing P R -f-)- Pr in the quark current. 
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A. B -> Kfj+fj,- 


In calculating the amplitude for the B —> Kfi + fi~ decay it is convenient to group the combinations of Wilson 
coefficients multiplying the same hadronic matrix element. Namely, the operators O i_6 mix at leading order into 
CV, 8,9 and it is customary to define effective Wilson coefficients as [37]: 


= 

CTo (Mb) = 


4tt 1 4 20 80 

- C 7-zC3-yC 4 --C 5 --C 6 , 

— C 9 + Y{q 2 ), 
a s 


47T 


a s 


C 


10 > 


^7,8,9,10(7*6) _ ^7,8,9,10 > 

a s 


( 5 ) 


where the function Y(q 2 ) at NLL can be found in Ref. [38]. We also incorporate the NNLL mixing of 0i and O 2 into 
O 7 and Og as calculated in Ref. [39]. The Wilson coefficients on the right-hand sides are evaluated at fi = fib- For the 
sake of readability we will from here on discuss only the effective Wilson coefficients that will be addressed simply as 
“Wilson coefficients” and denoted without the “eff” label. The values of the SM Wilson coefficients at scale fib are 
C 7 = -0.304, C 9 = 4.211, and C 10 = -4.103 [37, 38, 40], 

The decay spectrum as a function of the invariant mass of the muon pair is given by 

dT 2 

^2 (-B -> Kfi + fT) = 2 a^q 2 ) + -c^(q 2 ) , (6) 


where q 2 = +p M +) 2 , while functions a M (g 2 ), c AI (g 2 ) are combinations of Wilson coefficients and hadronic form 

factors and their explicit expressions can be found in Ref. [23] and in the Appendix of the present paper in the limit of 
m-n —>• 0. The rate depends on the sums of the Wilson coefficients of opposite chiralities, C 7 + C' 7 , Cg + Cg, Ci 0 + C [ 0 , 
from what follows that even in principle we cannot determine the chirality of the quark-current in B —> Kfi + fi~. 

Definitions of the hadronic form factors are relegated to the Appendix. We employ the form factors calculated in 
the unquenched lattice simulation using non-relativistic formulation of the b quark and staggered fermion formulation 
for the light quarks [24] . We use the z-expansion to parameterize the form factors and take into account the statistical 
errors given by the covariance matrix of the parameters, both given in [24]. However, we neglect additional systematic 
errors that should come on top of the ones contained in the covariance matrix. The correlations between form factor 
parameters are propagated onto observables of interest, namely we can construct y 2 statistic for B(B —> Kfi + fi~) 
and Rk , that are functions of the form factor parameters, as well as the Wilson coefficients. Nonlocal contributions 
to the decay amplitude due to operators 0 i ,2 are taken into account by leading order in operator product expansion 
together with next-to-next-to-leading logarithmic QCD corrections [39] . Higher orders in operator product expansion 
have been shown to have small effect in the large q 2 region [41]. Since the partial branching ratio that we are interested 
in corresponds to an integral over a large region of q 2 we rely on the semi-local quark-hadron duality [42]. In the SM 
limit the prediction of the branching ratio in the high-g 2 bin is 

B(B+ -a A'VV-)l* 2 M e[15 , 22 ]GeV2 = (10.2 ± 0.5) x 10- 8 . (7) 

The LHCb collaboration measured partial branching fractions below and above the region of charmonium reso¬ 
nances. For the q 2 > 15 GeV 2 region we can predict the partial branching ratio using form factors determined on the 
lattice that are largely free from extrapolation errors and parameterization dependence. Thus we will use [43], 

B(B+ K+fi+fi-) | 9 2 e[15j22]GeV 2 = (8.5 ± 0.3 ± 0.4) x 10” 8 , (8) 

as an experimental constraint, where the errors quoted are statistical and systematic, respectively. In our analysis we 
will sum the two and treat the observable with a Gaussian y 2 . 


B. B s —> fi + fi 

This decay receives contributions from operators with axial, scalar, and pseudoscalar lepton currents, and, owing 
to the pseudoscalar nature of the B s meson, the wrong-chirality Wilson coefficients will affect the decay with opposite 
sign. In the absence of (pseudo)scalar operators, the amplitude is proportional to the difference C 10 — C' w : 

P= 2 ^(C 10 -C[ 0 ), 

m Ba 


( 9 ) 
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and the “theoretical” branching ratio is expressed as 


B(B s ^ f i + fx-) th = Bo\P\ 2 , 


= fB s m B s G 2 F a 2 \V tb V t s\ 2 l~ 4m^ 
r s (4 tt)3 y m% t 


( 10 ) 


For the decay constant of the B s meson we take fs a = (228±8) MeV, consistent with the average made by FLAG [44]. 
Due to B s — B s oscillations and relatively large y s = Ar s /(2r s ) in the B s sector, the measured branching fraction 
actually corresponds to a time-integrated rate of the oscillating B s system to /i + n~ [45]. In effect, the value reported 
by the experimentalists is different from B(B S —> /t + /t _ ) th : 

B(B S -»■ M + M”) exp = T~~~2 [\ p \ 2 + 2/ S Re(P 2 )] . (11) 

1-2/s 

Latest average of the LHCb and CMS measurements of B s —> /t + /t _ branching fraction is [3] 

B(B S -»■ M + /V) exp = (2.8tg;5) x 10 ’ 9 ■ ( 12 ) 


The relative decay width difference y s = 0.061 ± 0.009 has been determined from LHCb simultaneous measurement of 
total width T s and width difference AT S in decay channels B s —> J/ipP + P~ [46]. The above determined value agrees 
very well with the HFAG and PDG averages [30, 47]. In the fits we use the values for P s and AT S reported by LHCb 
with summed statistical and systematic errors 

Ar s = (0.0805 ± 0.0123) ps" 1 , T s = (0.6603 ± 0.0042) ps" 1 , (13) 


with correlation coefficient —0.3 [46]. 


III. NEW PHYSICS IN C'g = -C \ 0 AND PREDICTION FOR R K 

We focus now on the SM extensions that affect the effective Hamiltonian solely by a single operator that is a 
product of right-handed quark current with a left-handed lepton current. In our operator basis it corresponds to a 
linear combination 0' 9 — 0' 19 implying 

C'(A) = -C( 0 (A), (14) 

where A is a scale where NP degrees of freedom are integrated out. An explicit example of such a scenario can be 
made in a leptoquark model that will be discussed in Section IV. If Eq. (14) holds at scale A it is neccessary to run the 
Wilson coefficients down to the low scale /if, using the renormalization group equations. Under QCD renormalization 
group the two operators do not run, keeping the constraint (14) intact [48]. 1 Thus we have, at low energies, a SM 
modification that satisfies 


C' 9 = -C[ o , (15) 

where C 9 10 are scale invariant, modulo small QED corrections. 

In Fig. 1 we show in gray the let region in the C' 10 plane as obtained from the fit to the partial branching fraction 
of B + —»• A' + /f + /r _ , cf. Eq. (8). The la region is defined here as y 2 < 2.30. The width of the “donut” reflects both 
experimental and form factor uncertainties. The SM point in the parameter space is marked with a dot and exhibits 
a tension with the measurement with \ 2 = 3.9. In Fig. 1 the lrr region (defined as before) of fit to the B(B S —► /U/t - ) 
according to Eq. (12) is depicted in blue. In this case the SM point is in comfortable agreement with the observable 
(\ 2 = 0.7). Then we perform combined fit to all of the above quantities and find the best value to be x 2 nin = 2.26, 
which is substantially better than the SM point with Xg M = 4.6. The green patch is defined by y 2 < Xmin +1 = 3.26 
with 39% C.L. and corresponds to the la region of predicted Rk given below. 

Assuming that the effective Hamiltonian (3), tailored for b —> se + e~, receives only SM contributions, unlike 
b —> sfi + that also receives NP contributions from C 9 10 , we can now predict the value of Rx- In Rk the 
uncertainties of the hadronic form factors cancel out to a large extent in the ratio and the formula boils down to: 

R k {C' w ) = 1.001(1) - 0.46 Re[C( 0 ] - 0.094(3) Im[C( 0 ] + 0.057(l)|C[ o | 2 . (16) 


1 Eq. (14) is broken only by tiny effects from QED renormalization. 
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Figure 1. Regions in the complex C( 0 plane that are in la agreement with B a — > q + fi (blue), B — > I\yA fj, (gray). Green 
area corresponds to the la coverage of Rk from fit to both observables. Black dot is the SM. 


Remaining uncertainties are indicated by the numbers in parentheses. In Fig. 2 we show contours of constant Rk 
in the C[ 0 plane using the formula (16) with central values for the coefficients. By dark gray we indicate the region 
corresponding to the measured value of Rk- In the same figure we plot again the 1 a prediction of C [ 0 , also shown 
in Fig. 1. We see an appreciable overlap with the measured Rk- Mapping the fitted region (green) to Rk we obtain 
the prediction 


R% ed - = 0.88 ± 0.08 , 

which is indeed in good agreement with 7?^- HCb = 0.745 ± 0:074 ±0.036 [1]. 


(17) 


A. Impact on B —> K*£ + £ 

B -> I<*l+£- is particularly interesting for the NP searches because of the observables that one can construct from 
the g 2 -dependent coefficients 1 1...9 which appear in the angular distribution, 

If sin 2 9k + If cos 2 9k + (If sin 2 9k ± If cos 2 9k) cos 2 9g 

+ I 3 sin 2 9k sin 2 9 c cos 2 (f> + I 4 sin 2 9k sin 29 c cos <f> + 15 sin 2 9k sin 9c cos 4 > 

+ (Ig sin 2 9k + If cos 2 9k) cos 9g + sin 2 9k sin 0 g sin </> (18) 

+Ig sin 29k sin 29c sin (f> + I 9 sin 2 9k sin 2 9c sin 2 cf) . 

Differential decay rate is then simply dr/ dq 2 = (37° + 67* — If — 27f )/4, and the similar expressions can be written for 
the transverse/longitudinal part of the decay rate, for the forward-backward asymmetry, A^q 2 ) = 3II/(4dT/dq 2 ), 
CP-asymmetry, and several other observables. Each of the coefficient functions, Ii = Ii(q 2 ), can be written in terms of 
transversity amplitudes, 7ljJ|| 0 (q 2 ), which are related to the respective spin states of the on-shell AT*-meson, and the 

amplitude A^' R (q 2 ) which is related to the off-shell virtual gauge boson decaying into the lepton pair. The superscripts 
L,R indicate the chirality of the lepton. Detailed expressions can be found, for example, in Refs. [38, 49-51]. 


d 4 T(B° -a K*°l+£~) _ 9 
dq 2 d cos 9g dcos9 K dcj) 32n 
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-0.5 0.0 0.5 1.0 

R e C 10 


Figure 2. Contours of constant Rk are indicated by dashed lines. Gray region represents the la measured range of Rk 
projected onto the C[ 0 plane, whereas green contour denotes the region allowed by B s —» n + and B —» K[i + fi~. Black dot 
is the SM. 


The strategy of looking for the NP effects through a detailed analysis of the angular distribution of B — > K*£ + £~ 
is somewhat plagued by hadronic uncertainties. The observables built up of A^^(q 2 ) turn out to be less sensitive to 
hadronic uncertainties because they involve the (combinations of) hadronic form factors which appear to be under a 
rather good theoretical control, especially in the region of small q 2 's [52-54] (see also discussion in Ref. [51]). On the 
other hand, the observables made of A u t '( q 2 ) entail the hadronic form factors that are less well understood. Moreover, 
the latter observables are subject to another kind of hadronic uncertainty, i.e. the one arising from misidentification of 
the Kn pairs coming from B —> K *(—>■ Ktt)£ + £~ with those emerging from B —> Ktt )£ + £~, where I\q stands 

for a broad scalar state [55-59]. Finally, and to avoid problems of the cc-resonances in the g 2 -spectrum of the decay, a 
standard strategy is to either work at low q 2 < m 2 ^, or large q 2 > 15GeV 2 , in which the impact of the cc-resonances 
is expected to be small. To be more specific, we fully rely on quark-hadron duality since we avoid the region in which 
the prominent narrow cc-resonances appear, and integrate over a window > 5 GeV 2 . 2 

With the information obtained in the previous section of this paper, i.e. with C[ 0 = —Cg extracted from the 
comparison of the measured B(B S —> fx + fx~) and B(B —> K^ + fi~) q 2 >15 Q e v 2 with the corresponding theoretical 
expressions, we already showed that we were able to verify the consistency of our result for Rk with the one measured 
at LHCb. With our approach, in which only the decay to muon-pair is modified, we can also predict Rk *, defined as 


T(i? —> K*fi + [i ) g 2 e [ 16 ] GeV 2 

r(-B —► /'f*e+e _ ) 9 2 6 [ 16 ] GeV 2 


(19) 


as well as the ratio of the two [6, 22], namely, 


Xk 


Rk* 
R k 


- 1 . 


( 20 ) 


In Ref. [19] it was shown that the ratio of forward-backward asymmetries integrated between q 2 £ [4,6] GeV 2 can 


2 For a recent attempt to more realistically model the effects of such resonances see Ref. [60] or those discussed previously in Refs. [61, 62]. 
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also be sensitive to lepton flavor universality violation. After defining, 

GeV 2 

I I 6 V) dq 2 

- " 


/4 GeV 2 


fbi 4 - 6 ] 4 r(s K*in-) q . em GeV2 ’ 

the ratio of forward-backward asymmetries is then simply, 


-Rfb — 


AV 

^^[4-6] 


A: 


fb[4—6] 


( 21 ) 


( 22 ) 


To compute the above-mentioned quantities we use the standard values of the Wilson coefficients [40], and include 
the effect of quark loops in the coefficients 67,9 arising from the operators Oi i2 , as calculated in Ref. [39]. We neglect 
the soft gluon corrections to the charm quark loop at low q 2 , which according to Ref. [63] is reasonable. At low q 2 the 
hard scattering contributions are neglected. For the form factors we use the values computed by means of QCD sum 
rules on the light cone [64], In Fig. 3 we show our results for Rk*, Xk and i?fb as functions of Re[C( 0 ]. For an easier 
comparison, in the same plot we also show Rk- The range 0.075 < RefC”^] < 0.41 has been obtained in the previous 
section of this paper, where we showed that for a given value of 0.075 < Re[C%] < 0.41 there is a region of allowed 
Im[Ci„], and therefore instead of curves in Fig. 3 we actually have the corresponding regions of values determined 
by Im[ 6 '] 0 ]. We should emphasize again that the uncertainties related to form factors cancel to a large extent in the 
ratios. As for the results, we first see that in the scenario with C' 10 = — Cg ^ 0, allowing coupling to muons only, and 





Figure 3. Rk, Rk* , Xk and Afb, defined in Eq. (1,19,20,22) respectively, are plotted as functions of RefCJo], in the range 
allowed by the measured values of B(B S — > fiAq~) and B(B —> K^ + fj,~) q 2 >15 GeV 2 - Instead of a curve for each quantity we 
actually have a region of values, reflecting the fact that for each Re[C(o] there is a range of allowed values of Im[Cio], as shown 
in Fig. 1. 

explicitly realized in the model with a (3, 2,1/6) leptoquark state, we get 


Rk = 0.88 ± 0.08 
X K = 0.27 ±0.19 


R k , = 1.11 ±0.08, 
i?fb = 0.84 ±0.12, 


( 23 ) 









which are obviously different from the values obtained in the SM, R^ 1 = 1.00, = 0.996(5) ss 1, i?fb = 0.995(4) « 

1, and Xk = —0.004(5) ss 0. Notice, however, that while our value for Rk is lower than the one in the SM, our 
prediction for Rk* is larger than that obtained in the SM. The measurement of Rk* at LHCb will therefore help to 
either confirm or discard our model as a viable description of the lepton flavor universality violation. The errors in 
Eq. (23) are completely dominated by the range of RefC^g] and Im[G" 0 ], while those arising from form factors are 
reduced in the ratios and induce an uncertainty negligible in comparison with that coming from the variation of C [ 0 . 

Besides the above quantities, one can also check on the asymmetries A^ and Ay m \ defined in the Appendix, 
which are experimentally more difficult to study but which could be very useful to compare with predictions as their 
values can considerably change if Cg 10 ^ 0. To exemplify that feature we consider the bin q 2 £ [2, 6] GeV 2 and in the 

SM we have (A^^g 2 ))™ j 2 Gey2 = —0.05(2), and (Ay n b(9 2 ))p4 g j Gey2 ~ 0, for either electrons or muons in the 
final state. If, instead, the coefficients C[ 0 = — Cg become non-zero and take values within the green region shown in 
Fig. 1, then in the case of B —» I\*n + , the above values change to 


<4 2 V)> 


g 2 e[2,6] GeV 2 


~ -0.09 


,-0.24 


for Re[C( o ]=0.08 


| (A 


for Re[C( o ]=0.4 


(Im) 

T 


(Q ))q2 e [2,6] GeV 2 I 


< 0.27. (24) 


Notice also that (Ay 111 ^ 2 )) has not yet been measured, and that the current errors on {A^\q 2 )) are still too large 
for making a meaningful quantitative comparison with our results [65]. 

Finally, before closing this part of our paper, we need to comment on P/(g 2 ), an observable constructed from 
coefficients of the angular distribution of the B —► K*( + £~ decay [66], P/(g 2 ) = Is/yJ—A /f/|, which has been 
measured at LHCb, and turned out to be 4 a away from the value predicted in the SM when integrated over an 
interval q 2 £ [4.3, 8.68] GeV 2 [67]. More specifically, the SM value is (Pg)^_ g g8 , = —0.90(5), while the measured one 

is ( P 5)f 4 H 3 - b 8.68] = -0.19(16) [67], which can be compactly written as, = 0.22(18). While 

the interpretation of this discrepancy is somewhat controversial [18, 19, 68], it is nevertheless interesting to check 
whether or not the leptoquark model used in this paper (and discussed in more details in the following Section) can 
describe the manifest disagreement between theory and experiment. With the values of C' w = —Cg discussed above 
we indeed see that (Pg^i-s 68] < 1> but with the leptoquark model discussed here we cannot reach 

very low values. We instead obtain 0.78 < (-Pg)j/3_8 68]/^ P s)f4 3-8 68] — 0-98. A similar tendency is observed for other 
bins, and in particular the one corresponding to q 2 £ [1,6] GeV 2 . 


IV. MODEL WITH A SCALAR LEPTOQUARK 

In this Section we discuss a specific model in which the scenario discussed above, i.e. Cg = — C' w , is explicitly 
realized and involves the presence of a light scalar leptoquark state A. More specifically, we choose the leptoquark 
A to carry the quantum numbers (3,2,1/6) of the SM gauge group. Its couplings to fermions are described by a 
renormalizable Lagrangian 

C = YijLi ir 2 A*dnj + h.c. 

= Y zj (-l Li d Rj A< 2 / 3 >* + v Lk (V™™)ld Rj M-^*) + h.c., (25) 

where Y is a 3 x 3 complex matrix, Lj and d R) are the lepton doublet and down-quark singlet. Charge eigenstates 
of the leptoquark doublet are denoted with A^ 2 / 3 ) and A* -1 / 3 ) and we will assume that they are degenerate. The 
second line in the above Lagrangian is written in the fermion mass basis, and a relative PMNS rotation in lepton 
doublet components has been assigned to the neutrino sector. 

Clearly, the lepton flavor universality is explicitly broken by the terms presented in Eq. (25). This might appear 
questionable because in a similar situation in which the coupling of leptoquark to ^ic would be allowed, the ratio of 
the electronic and muonic widths of the decay of J/'if) and its radial excitations have been accurately measured, and 
shows no violation of the lepton flavor coupling universality. In particular, the measured Y(J/ip —> —> 

e + e~) = 1.0016 ± 0.0031 [30] is in excellent agreement with its SM value, 1.00001. 3 That situation is, however, 


3 By explicitly including the lepton mass in the calculation of phase space we obtain, 

167ro 2 ( 2m? 


r (J/ip -s. £+£-) = 


27 m 


j/i/j 


1 + 




1 - 


4m 2 

m J/* 


and the effect on the ratio of the electronic and muonic widths is extremely small. 


fj/ tp 
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much different from the examples discussed in this paper, because the amplitude for J/ip —> £ + £~ is dominated by 
the tree-level electromagnetic interaction diagram which is much larger than the weak interaction one, suppressed by 
l/m| with respect to the dominant one, and therefore completely negligible. Our leptoquark state is TOa mz, and 
its contribution to J/ip —» £ + £~ is even smaller than the weak interaction diagram and cannot make an impact on 
the decay of charmonia at the present level of accuracy. 

Instead, the weak b —> s/j, + fi~ decays in the SM are loop-induced so that the tree level contribution involving 
couplings to the leptoquark state may become comparable in size to the SM amplitude, which is why the b — > s/i + /i _ 
is likely to be more sensitive to the presence of the term described by the lagrangian (25). The relevant leptoquark 
coupling for the b —► s/t + /t“ is the product Y^Y^g, which enters the Wilson coefficients divided by m 2 A . The scalar 
particle exchange generates scalar operators in the Fierzed basis and those appear as (pseudo)vector currents in the 
ordinary operator basis [9]: 


Cl ° ° 9 2V2G F V tb V t * s a 


Y uY * 

1 /ib 1 


m 


A 


(26) 


We assume other elements of Yukawa matrix Y to vanish. The same state will also contribute at loop level to electro- 
and chromo-magnetic operators C' 7 {m&) and C'g(mA) where these coefficients will be suppressed by electromagnetic 
a(wA)/(47r) and strong as(mA)/(47r) couplings at high scale ttia, respectively. We have explicitly checked that these 
modifications result in negligibly small value of C 7 when compared to the C 7 of SM, cf. Eq. (5). In the remainder of 
this Section we will analyze additional observables that constrain this leptoquark scenario. 

The considered leptoquark state A couples to the neutrinos with the same couplings as to the charged leptons, only 
modified by a PMNS rotation matrix. Namely, the charge —1/3 state will generate ( sb){vv ) operators while the box 
diagrams will lead to B s — B s mixing. 


A. Contribution of (3,2,1/6) leptoquark in B s — B s oscillation frequency 


The state (3, 2,1/6) will induce A B = 2 box diagrams with /i and A^ 2 / 3 ) or v and A^ running in the box. The 
two contributions of boxes with [i and v are equal in the m^ = 0 limit and in sum they amount to 


C 6 LQ (m A ) = 


1 fib 1 /is 


647r 2 m^ 


(27) 


The effective A B = 2 Hamiltonian is defined as 


n cS = Cf M {b lfi P L s) (b^P L s) + C^ Q (h„PRs) (b^PRs) + h.c., 


(28) 


where Pl/r = (1 T 7s)/ 2- The coefficient in Eq. (27) is subject to QCD renormalization and has to be evaluated 
at scale The anomalous dimensions of Cg® is however equal to the one of Cf M . Therefore the two Wilson 
coefficients renormalize with the same multiplicative factor between scales /r = m t , where SM is matched onto 
effective Hamiltonian (28), and /if,, where the hadronic matrix elements are computed. Remaining Cg ^ running from 
toa down to m t is already in the asymptotic regime of QCD and can be safely neglected. The mass difference of the 
B s — B s system is then 


A mB g = 



Gpm^ 

16ir 2 


(V; b V ts ) 2 VB S 0 (x t ) + 


4 


Cg Q (TOA) 


-75)s&7 /i (l -75 )s|£°). 


(29) 


By using Eq. (26) we can write 


C 6 LQ (m A ) = -^(VtlV t g) 2 a 2 m 2 A (C [* 0 ) 2 , 

which, together with - -/ 5 )sb"/^(l - 7 s)s|B°) = (8/3)/| m| B Bs , gives 


Ator. = 


r f ,,l w 

67T 2 


\V t * b V t g\ 2 f 2 m B3 B B ,T) B So(x t ) 


1 - 


1 


Q 2 rr /^2 TO A 
27T 2 S'o(xt) 10 rn^Y 


Am™ 


(30) 


(31) 
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With the current values for fg, = 228(8) MeV and Bb s = 1.33(6), as obtained in numerical simulations of QCD on 
the lattice [44], and m^ s (?n t ) = 160^ GeV [30], we get 4 

Am® = 17.3 ± 1.7 ps -1 , (32) 

which is in excellent agreement with the measured A ms s = 17.7(2)ps _1 [30]. With the values of C' w determined in the 
previous Section, we see that Eq. (31) leads to a very loose upper bound for toa- For example, for Re[C( 0 ] £ [0.15, 0.35], 
we get the upper bound of the order 100 TeV. 


B. Impact of (3,2,1/6) leptoquark on B —> Kvv 

In the presence of leptoquark A the pair of neutrinos in the final state of B —> Kvv may be in any flavor combination. 
In order to encompass such a possibility we must extend the effective Hamiltonian of Ref. [70] to account for the 
disparity in neutrino flavors: 


AC 

n cS = + c%o%) . (33) 

The operators are defined as O l [ R = (sj M PL,_Rb)(vij M (l — 75 )vj). The authors of [70] found that in the SM the 
Wilson coefficient at next-to-leading order in QCD is 

Cf M = Cf = -6.38 ± 0.06 , (no sum over i implied). (34) 


If the leptoquark state (3, 2,1/6) is present then it will manifest itself in B —> Kvv through right-handed operators: 


c% = 


1 (VY) ib (VY)* ja 
N 4m^ 


N _ G F V tb V* s a 
spin 


(35) 


Here V denotes the PMNS matrix. The experimentally accessible decay width of B —>• Kvv is a sum of partial widths 
of B —^ KCpi. The amplitudes are proportional to the sum of the SM and leptoquark contribution and the two will 
interfere in the B ^ Kvv decays width as 


T(B Kvv) ~ ^ | S XJ Cl M + 
i,j =1 

= 3|Cf M | 2 -f |C ( 0 | 2 — 2Re[C'f M *C'( 0 ]. 


(36) 


C[ 0 is the Wilson coefficient of b —> that we obtained from the fits to experimental data in the previous 

Section. Last line of Eq. (36) was obtained by applying the unitarity of matrix V, and assuming that Y^b and Y^ s 
are the only non-zero elements of the matrix Y. Finally, the g 2 -spectrum of this decay reads, 

^(B -a Kvv) = 3 g ^j^ f+(q 2 ) [\(m 2 B ,m 2 K ,q 2 )f 2 (^3|Cf M | 2 + |C ( 0 | 2 - 2Re[G| M *C( 0 ]) , (37) 

where q 2 in this case stands for the invariant mass of the neutrino pair. Notice that the above expression, for C' 10 = 0, 
confirms Eq. (2.14) of Ref. [70]. The expression (37) can be recast into a product of the SM < 7 2 -spectrunr and a 
correction factor, 


1.01 < 


'4 


c' w /c\ 


SM 


I 2 - l MC[o/C : l 


SMi 


< 1.05 


(38) 


where its lower and upper bounds have been derived from the ler region of C( 0 , obtained in the previous Section. We 
learn that the B(B —> Kvv) will increase by at most 5% if leptoquark A is present. 


4 To evaluate Ara® we also used rjg = 0.55(1) [69], and So(xt) = 2.25/],/, the Inami-Lim function at xt = 
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V. SUMMARY AND CONCLUSION 

In this paper we discussed a possibility of constraining a scenario of New Physics affecting the b -A sp + p,~ decays 
through coupling with the operators O 9 10 . Such a scenario is explicitly verified in a model with a light scalar 
leptoquark state, A, carrying the quantum numbers (3,2,1/6) of the Standard Model gauge group. In this scenario, 
Cg = —C [g is specific for the muons in the final state. In the leptoquark model discussed in this paper, C' w is related 
to Y^Y* s /m\. From the currently available experimental data on B(B S -a p + p~) and B{B -a Kp + /r _ )[i 5 i 22 ]GeV 2 
we were able to constrain Re[Cj 0 ] and Im[C'( 0 ], which are then used to compute Rk = B(B -A Kp, + p~)/B(B -a 
Ke + e~)[i t 6 jQ eV 2 = 0.88(8). We find a good agreement with a recent experimental result, i?^- HCb = 0.75(14). 

After having passed this test we focused on B A K*t + £~ decays. Within our scenario, and the range of C' w 
obtained from B(B S -A p + pT) and B(B -A Kp + pL~)y i 5i 22 ]GeV 2 > we predict Rk * = B(B -A K*p + p~)/B(B —t 
A'*e + e _ )[ 16 ] GeV 2 = 1 . 11 ( 8 ), and Rk*/Rk = 1.27(10). Similarly, in this scenario the ratio of forward-backward 
asymmetries becomes different from unity. In particular, we find Rfb = GeV 2 = 0.8(1). Furthermore, 

we checked that a combination of coefficients of the angular distribution of B —> K* p, + p~, known as P 5 , and weighed 
over a specific bin of q 2 ’ s, indeed becomes smaller than its value predicted in the Standard Model. However, it cannot 
explain a very low value of |(Ps)|[ 4 . 3 - 8 . 68 ] measured at LHCb, for which the Standard Model prediction is still a 
subject to controversies mainly related to the issue of treatment of the charm quark loops. 

Finally, in the leptoquark model our constraints on the Wilson coefficient C [ 0 can have impact on other physical 
processes. We checked, in particular, that the contribution to the frequency of oscillation in the B s — B s system is 
insignificant, and that only up to five percent enhancement of B{B -A Kuv) can be obtained. 
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Appendix A: B —> form factors 


For completeness we remind the reader of the standard parameterization of the B 
elements in terms of the relevant form factors, 


hadronic matrix 


{K{k)\sj^b\B(p)} 

(K(k)\scr^b\B(p)) 


ip + % 
-i {p^K - 


Tor — m 


K . 




/ + (<? 2 ) 


Pvk^) 


2/r(g 2 ) 
ms + rriK 





<3VJo(<7 2 ), 


(Al) 


(K*(k)\s'y fl (l - ^ 5 )b\B{p)) = e M „p 0 -e*V’fc' 

A 2 (q^ 


2 V(q 2 


+i{p+k)^(£* ■ q) 


ms + rriK * 

2 rriK* 


ms + rriK * 


+ iq^(e* ■ q) 


ie*(m B + rriK*)Ai{q 2 ) 
[A 3 (q 2 ) — Ao(g 2 )] , 


(A2) 


{K*(k)\sa^ v q v {l +j 5 )b\B(p)) 


^£/ii2pcr£ 

+ [<(™ 

+ (£*W) 


* v pPk a Ti(q 2 ) 

b - ™ 2 k*) - (e* ' Q) (2 P~q)p] T 2 (q 2 ) 
0,2 

„2\ 


% 


m B — m 


(P + k )p 


K* 


W 


(A3) 


where 2 niK*A 3 (q 2 ) = (ms + nriK*)Ai(q 2 ) — (ms — rriK*)A 2 (q 2 ), and Tj(0) = 72(0). In the limit of massless lepton, 
the g 2 -dependent functions entering eqs. (6,18), relevant for the present study, read 


= IS = \ (l^’T + l^’T) , 

I c 2 = ~\A L 0 ' R \ 2 , I s 6 = 2 Re (A b Af * - A«A* *) , 


(A4) 
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where 


n<?) 


NX(q 2 ) 


(Cg + C 9 )f + (q )+ (C7 + C 7 )/t(<7 ) 

2 

+ 

2" 

(C 10 + C'( 0 )/ + (<? 2 ) 

m B + m K 




^’V) 

A\( R {q 2 ) 

A^ R (q 2 ) 


\J Nq 2 X(q 2 

-V N q 2 (m 2 B - m 2 K .) 
sfN 


^(C 7 + C' 7 ) Ti(q 2 ) + [(C 9 + Q =F (C w + C[ 0 )\ 


V(q 2 ) 


ms + mx * 


^?(CV - C' 7 ) T 2 (q 2 ) + [(C 9 - C' 9 ) T (C w - Cio)] 


Al(q 2 ) 

m B - m K * 


2niK * 


[(C 9 - C' 9 ) =F (C 10 - Cio)] X 


(m B - m\, - q 2 )(m B + m K *) A 1 (q 2 ) - A (q 2 ) 


A 2 (q 2 


m B + m.K * 


+ 2 m& (C7 — Cf) 


(m 2 B + 3 m 2 K , - q 2 ) T 2 (q 2 ) - 


X (q 2 ) 


~T 3 (q 2 ) 




(AS) 


and N = \V t bV t * s \ 2 X 1 11 ^ 2 (q 2 ) Gpa 2 /(3072ir 5 m B ), X(q 2 ) = [q 2 — ( m B + tuk*) 2 ] [q 2 — ( m B — ttik*) 2 ]- In the same 
massless lepton limit, a p = —c M , If = 3/f and If = —If, so that dr(i? —> Kfj,fi)/dq 2 = 4a M /3, and dT(B —> 
K* fin)/dq 2 = 4 If —If = \A ^ R \ 2 + \A ^' R \ 2 + \Aq' R \ 2 . The two transverse asymmetries discussed in the text are defined 

as 4 2 V) = J 3 /(2 If) = (|AJ 2 - \A^ 2 )/{\A 1 _\ 2 + lAyl 2 ), and 4 Im) (g 2 ) = I 9 /(2If) = -2Im(A || Al)/(|A ± | 2 + |An| 2 ). 
In the computation of decay amplitudes we rely on the full QCD form factors computed on the lattice (in case of 
B — ► Kit) or by means of the QCD sum rules near the light cone (in case of B —> K*M) and to the standard (universal) 
Operator Product Expansion. We did not rely on the effective theory approaches such as QCD factorization because 
we do not know how to reliably compute the relevant form factors and include all power corrections. At the level of 
present accuracy however, both approaches lead to compatible results. 


[1] R. Aaij et al. (LHCb), Phys.Rev.Lett. 113, 151601 (2014), 1406.6482. 

[2] G. Hiller and F. Kruger, Phys.Rev. D69, 074020 (2004), hep-ph/0310219. 

[3] V. Khachatryan et al. (CMS, LHCb), Nature 522, 68 (2015), 1411.4413. 

[4] C. Bobeth, M. Gorbahn, T. Hermann, M. Misiak, E. Stamou, et al., Phys.Rev.Lett. 112, 101801 (2014), 1311.0903. 

[5] R. Alonso, B. Grinstein, and J. Martin Camalich, Phys.Rev.Lett. 113, 241802 (2014), 1407.7044. 

[6] G. Hiller and M. Schmaltz, Phys.Rev. D90, 054014 (2014), 1408.1627. 

[7] S. L. Glashow, D. Guadagnoli, and K. Lane, Phys.Rev.Lett. 114, 091801 (2015), 1411.0565. 

[8] B. Bhattacharya, A. Datta, D. London, and S. Shivashankara, Phys.Lett. B742, 370 (2015), 1412.7164. 

[9] N. Kosnik, Phys.Rev. D86, 055004 (2012), 1206.2970. 

[10] D. Ghosh, M. Nardecchia, and S. Renner, JHEP 1412, 131 (2014), 1408.4097. 

[11] S. Biswas, D. Chowdhury, S. Han, and S. J. Lee, JHEP 1502, 142 (2015), 1409.0882. 

[12] T. Hurth, F. Mahmoudi, and S. Neshatpour, JHEP 1412, 053 (2014), 1410.4545. 

[13] W. Altmannshofer and D. M. Straub (2014), 1411.3161. 

[14] S. Sahoo and R. Mohanta (2015), 1501.05193. 

[15] I. d. M. Varzielas and G. Hiller (2015), 1503.01084. 

[16] B. Gripaios, M. Nardecchia, and S. Renner (2014), 1412.1791. 

[17] S. Jager and J. Martin Camalich, JHEP 1305, 043 (2013), 1212.2263. 

[18] S. Descotes-Genon, J. Matias, and J. Virto, Phys.Rev. D88, 074002 (2013), 1307.5683. 

[19] W. Altmannshofer and D. M. Straub, Eur.Phys.J. C73, 2646 (2013), 1308.1501. 

[20] A. Crivellin, G. D’Ambrosio, and J. Heeck (2015), 1501.00993. 

[21] S. D. Aristizabal, F. Staub, and A. Vicente (2015), 1503.06077. 

[22] G. Hiller and M. Schmaltz, JHEP 1502, 055 (2015), 1411.4773. 

[23] D. Becirevic, N. Kosnik, F. Mescia, and E. Schneider, Phys.Rev. D86, 034034 (2012), 1205.5811. 

[24] C. Bouchard, G. P. Lepage, C. Monahan, H. Na, and J. Shigemitsu (HPQCD Collaboration), Phys.Rev. D88, 054509 
(2013), 1306.2384. 

[25] R. R. Horgan, Z. Liu, S. Meinel, and M. Wingate, Phys.Rev. D89, 094501 (2014), 1310.3722. 

[26] I. Dorsner, S. Fajfer, and A. Greljo, JHEP 1410, 154 (2014), 1406.4831. 

[27] G. A ad et al. (ATLAS), Eur.Phys.J. C72, 2151 (2012), 1203.3172. 

[28] Tech. Rep. CMS-PAS-EXO-12-042, CERN, Geneva (2013), URL http://cds.cern.ch/record/1542374. 

























13 


[29] M. I. Gresham, I.-W. Kim, S. Tulin, and K. M. Zurek, Phys.Rev. D86, 034029 (2012), 1203.1320. 

[30] K. Olive et al. (Particle Data Group), Chin.Phys. C38, 090001 (2014). 

[31] S. Weinberg, Phys.Rev. D22, 1694 (1980). 

[32] W. Buchmuller and D. Wyler, Pliys.Lett. B177, 377 (1986). 

[33] I. Dorsner, S. Fajfer, and N. Kosnik, Phys.Rev. D86, 015013 (2012), 1204.0674. 

[34] B. Grinstein, M. J. Savage, and M. B. Wise, Nucl.Phys. B319, 271 (1989). 

[35] M. Misiak, Nucl.Phys. B393, 23 (1993). 

[36] A. J. Buras and M. Munz, Phys.Rev. D52, 186 (1995), hep-ph/9501281. 

[37] A. Buras, M. Misiak, M. Munz, and S. Pokorski, Nucl.Phys. B424, 374 (1994), hep-ph/9311345. 

[38] W. Altmannshofer, P. Ball, A. Bharucha, A. J. Buras, D. M. Straub, et al., JHEP 0901, 019 (2009), 0811.1214. 

[39] C. Greub, V. Pilipp, and C. Schupbach, JHEP 0812, 040 (2008), 0810.4077. 

[40] C. Bobeth, M. Misiak, and J. Urban, Nucl.Phys. B574, 291 (2000), liep-ph/9910220. 

[41] B. Grinstein and D. Pirjol, Phys.Rev. D70, 1]L4005 (2004), hep-ph/0404250. 

[42] M. Beylich, G. Buchalla, and T. Feldmann, Eur.Phys.J. C71, 1635 (2011), 1101.5118. 

[43] R. Aaij et al. (LHCb), JHEP 1406, 133 (2014), 1403.8044. 

[44] S. Aoki, Y. Aoki, C. Bernard, T. Blum, G. Colangelo, et al., Eur.Phys.J. C74, 2890 (2014), 1310.8555. 

[45] K. De Bruyn, R. Fleischer, R. Knegjens, P. Koppenburg, M. Merk, et al., Phys.Rev. D86, 014027 (2012), 1204.1735. 

[46] R. Aaij et al. (LHCb), Phys.Rev.Lett. 114, 041801 (2015), 1411.3104. 

[47] Y. Amhis et al. (Heavy Flavor Averaging Group (HFAG)) (2014), 1412.7515. 

[48] C. Bobeth, P. Gambino, M. Gorbahn, and U. Haisch, JHEP 0404, 071 (2004), hep-ph/0312090. 

[49] C. Bobeth, G. Hiller, and G. Piranishvili, JHEP 0807, 106 (2008), 0805.2525. 

[50] S. Descotes-Genon, T. Hurth, J. Matias, and J. Virto, JHEP 1305, 137 (2013), 1303.5794. 

[51] D. Becirevic and E. Schneider, Nucl.Phys. B854, 321 (2012), 1106.3283. 

[52] J. Charles, A. Le Yaouanc, L. Oliver, O. Pene, and J. Raynal, Phys.Rev. D60, 014001 (1999), liep-ph/9812358. 

[53] M. Beneke, T. Feldmann, and D. Seidel, Nucl.Phys. B612, 25 (2001), hep-ph/0106067. 

[54] M. Beneke, T. Feldmann, and D. Seidel, Eur.Phys.J. C41, 173 (2005), hep-ph/0412400. 

[55] D. Becirevic and A. Tayduganov, Nucl.Phys. B868, 368 (2013), 1207.4004. 

[56] D. Das, G. Hiller, M. Jung, and A. Shires, JHEP 1409, 109 (2014), 1406.6681. 

[57] T. Blake, U. Egede, and A. Shires, JHEP 1303, 027 (2013), 1210.5279. 

[58] J. Matias, Phys.Rev. D86, 094024 (2012), 1209.1525. 

[59] M. Doring, U.-G. Meiner, and W. Wang, JHEP 1310, 011 (2013), 1307.0947. 

[60] J. Lyon and R. Zwicky (2014), 1406.0566. 

[61] A. Ali, T. Mannel, and T. Morozumi, Pliys.Lett. B273, 505 (1991). 

[62] F. Kruger and L. Sehgal, Phys.Rev. D55, 2799 (1997), hep-ph/9608361. 

[63] A. Khodjamirian, T. Mannel, A. Pivovarov, and Y.-M. Wang, JHEP 1009, 089 (2010), 1006.4945. 

[64] P. Ball and R. Zwicky, Phys.Rev. D71, 014029 (2005), hep-ph/0412079. 

[65] R. Aaij et al. (LHCb), JHEP 1308, 131 (2013), 1304.6325. 

[66] S. Descotes-Genon, J. Matias, M. Ramon, and J. Virto, JHEP 1301, 048 (2013), 1207.2753. 

[67] R. Aaij et al. (LHCb), Phys.Rev.Lett. Ill, 191801 (2013), 1308.1707. 

[68] S. Jager and J. Martin Camalich (2014), 1412.3183. 

[69] A. J. Buras, M. Jamin, and P. H. Weisz, Nucl.Phys. B347, 491 (1990). 

[70] W. Altmannshofer, A. J. Buras, D. M. Straub, and M. Wick, JHEP 0904, 022 (2009), 0902.0160. 



